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We re-examine recent determinations of power corrections from r-decay and confront the results with the existing ones from 
QCD spectral sum rules (QSSR). We conclude that contrary to the QSSR analysis, which lead to (a a G 2 ) = (6.8 ± 1.3)10~ 2 
GeV 4 , r-decay is not a good place for extracting the gluon condensate due to its extra a 2 coefficient which suppresses its 
contribution in this process. Results from e + e~ sum rules and r-decay: pa s {uu) 2 = (4.5 ±0.3) 10 -4 GeV 6 , where p = 3.0 ±0.2 
confirm the deviation from the vacuum saturation estimate of the four-quark condensate. "Non-standard" power corrections 
(direct instantons, duality violation and tachyonic gluon mass) beyond the SVZ-expansion, partially cancel out in the V±A 
hadronic r-decay channel, which gives at order a 4 : a s (M T ) — 0.3249 (29) cx (75)th leading to a s {Mz)\r = 0.1192 (4) cx (9)th, in 
remarkable agreement with (but more accurate than) a s (Mz)\z = 0.1191 (27) obtained at the same a 4 order from the Z-width 
and the global fit of electroweak data. Finally, the role of the tachyonic gluon mass in the determinations of \V US \ from r-decay 
and of fh a from r-decay, e + e~ and (pseudo)scalar channels is emphasized. 



1. Introduction 

Hadronic V + A r-decay is expected to provide the most 
accurate determination of a s once its value obtained 
at the r-mass is runned until the Z-mass |1I2I3I4I5] , 
The extraction of a S) at this relatively low scale M T , 
can become feasible, as theoretically, the V + A chan- 
nel is less sensitive to the QCD non-perturbative ef- 
fects than the individual V and A channel (if one uses 
the OPE a la SVZ [6]) and to the values of the quark 
and gluon condensates obtained from an overall fit of 
the hadronic channels [60, e+e" |8|9|10|11|12] and r- 
decay |13|14|15|16|17|18|19j data, light baryon sum rules 
[217] , heavy quarkonia sum rules |6I21I22I2"5] and heavy- 
quarkonia mass-splittings [23]. In addition, the domi- 
nant theoretical uncertainties on the a s determination, 
though relatively modest at the r-mass (10%), becomes 
tiny at the Z-pole as it decreases as 1/ log 2 . More quan- 
titatively, the non-strange AS — component of the 
r-hadronic width can be expressed as: 

_ r(r -> v T + hadrons| A s=o) 
T(r — > I + Di ± u T ) 

= i\V u d\ 2 SEW x 

( 1 + «S<°> + S' EW + 8$ + S svz + S ust ) , ( 1 ) 

where 5 m is the light quark mass corrections, 

8 

<5 SVZ = S{D) ' ( 2 ) 

is the sum of the non-perturbative (NP) contributions 
of dimension D within the SVZ expansion [6] , while <5 ns t 
are some eventual NP effects not included into S svz . We 
shall use: \V ud \ = 0.97418 ± 0.00027 gS]. The elec- 
troweak corrections are: 

Sew = 1-0198 ± 0.0006 [26] and S' EW = 0.001 HJ|.(3) 
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<5(°) is the perturbative correction, while S m is the light 
quark masses corrections. They will be discussed in de- 
tails in the next paragraph. 

2. SVZ power corrections in hadronic r-decay 

Using the QCD expressions compiled in BNP [2], the 
light quark masses corrections read: 




where fh q and a s = (a s /ir) are respectively the light 
quark running masses and QCD running coupling eval- 
uated at M T . Using the recent determinations of the 
running mass in units of MeV and evaluated at 2 GeV 
|28l29l30l3lj : 

fh s = 96.1 ± 4.8 , m d = 5.1±4, m u = 2.8 ± 0.2 , (5) 

one can deduce the small and negligible corrections in 
10~ 4 : 

€]v = "(3.15 ±0.51) , S% v+A = -(4.22 ±0.68) .(6) 

A theoretical estimate of the quark masses contributions 
has been often used as input in existing analysis of r- 
decay data [14115116] . The non-perturbative corrections 
a la SVZ [B] read to leading order in m q and a s [20: 

±167r 2 (m n =F m^^uipu ± $ d 4>d) 
« M^,-^ 2 ( as G 2 ) 3 (7) 

2 We shall use the complete QCD expression in BNP [2] for the 
numerics. 
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Table 1 

SVZ power corrections from r-decay compared with the ones from 
e~*~e — and QSSR analysis of the other hadronic channels. FO and CI 
correspond to a s (M T ) = 0.331 ± 0.013 and 0.350 ± 0.010. ALEPH 
fits come from [T4 16 , while OPAL fits arc from |15| . We take the 
average of different results and take the quadratic mean of the error 
(bold face) when the different fits are in good agreement, while in the 
case where some of the results are not significant, wc only consider 
the most accurate fit (boldface). The final sum <5 avz comes from the 
average or/and from the most accurate determinations. 
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where (jpip) = (uu) ~ (dd) and p ~ 2 - 3 |7l8l9ll0ll71 
118119120] indicates a deviation from the vacuum satu- 
ration of the four-quark condensate estimate. Vacuum 
saturation has been assumed for estimating the dimen- 
sion D=8 condensates which is expected to be a very 
rough approximation [8117118111)] . We compare in Ta- 
ble [TJ the size of the given dimension non-perturbative 
(NP) contributions obtained from a fit of the r-decay 
data with the ones obtained from QCD spectral sum 
rules (QSSR) [6 7J using the 1=1 component of the e + e~ 
data [8]. Some comments follow: 
• General remarks 

The results from QSSR have been obtained using the 
ratio of the Laplace/Borel (LSR)sum rules where radia- 
tive corrections tend to cancel out. Therefore, contrary 
to r-decay, the result will be only slightly affected by the 
truncation of the PT series. One can inspect in Table 
[T] that there is a fair agreement (within the errors) be- 
tween the different results from LSR used in e + e~ and 
from T-decay data analysis. One can notice that r-decay 



analysis can lead to inaccurate results for some NP ef- 
fects in the channels where they are tiny. However, the 
smallness of these NP effects in the V+A channel makes 
it a best place for extracting a s . 
• Gluon condensate 

An attempt to extract this quantity from r-decay in 
different channels |14|15|16] leads to unreliable results 
because its value which is expected to be universal dif- 
fers from one channel to another. In units of 10~ 2 GeV 4 , 
r-decay gives in the V, A and V+A channels: 

(a s G 2 )~ -(0.8 ± 0.4V, -(2.2±0.4) A , 

-(1.5±0.3V +A . (8) 

One can also worry on the negative sign which differs 
from the one obtained from different QSSR analysis 
[6I7I8I9I10I11I12I13I21I22I2"3 24fJ. This can be due to the 
fact that its contribution in these r-decay channels is dif- 
ficult to extract from the data as it acquires an extra a 2 
correction and then becomes relatively tiny compared 
to the other condensate effects. Indeed, from Eq. ([7]), 
one can notice that it is one order of magnitude smaller 
than the mfyip) quark and of the four-quark condensate 
effects. Another problem may arise that in ALEPH and 
OPAL analysis, one has to do simultaneously a fit of 
many parameters in the r-decay analysis (a s , D = 4, 6 
and 8 condensates). This is not the case of the QSSR 
analysis of e + e _ within LSR, which has a stronger sen- 
sitivity to (a s G 2 ) and then permits its robust estimate. 
In these QSSR analysis, it has been found from differ- 
ent analysis a strong correlation between the D = 4 and 
D = 6 condensate contributions. Using the average of 
the ratio of the gluon and four-quark condensates [8] 
determined in [9I10I11I12I13] : 

r 46 = {a ^\ = (106 ± 12) GeV" 2 , (9) 
pa s {uu) z 

which reduces the analysis to a one-parameter fit Q, the 
LSR analysis of e + e~ data gives [8]: 

(a s G 2 )| e+e - = (6.1 ± 0.7)10~ 2 GeV 4 , (10) 

while the heavy quarkonia mass-splittings [24] lead to: 

(a s G 2 )| hcavy = (7.5 ± 2.5)10- 2 GeV 4 , (11) 

from which we deduce the average: 

(a s G 2 ) = (6.8 ± 1.3)10~ 2 GeV 4 , (12) 

where the errors have been averaged quadratically. 
Though higher by about a factor 1.8 than the original 
SVZ value : 

(a s G 2 )| svz ~3.8xl0- 2 GeV 4 , (13) 

the previous estimate is at the border of the lower bound 
allowed by Bell & Bertlmann [21] in the analysis of the 

3 Using arguments based on magnetic confinement, one can argue 
that (q s G 2 ) is positive |32| . For a review, see e.g. I33| . 
4 The fit has been performed using a standard Mathematica least- 
square fit program. 
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heavy quark sum rules and is in line with results from 
Finite Energy Sum Rules (FESR) analysis [TO] or some 
other methods [11112] in e + e _ data and with the avail- 
able lattice calculations [34] 0: 



(a s G 2 )| 



latt 



(12.6±3.1)10~ 2 GeV 4 



(14) 



Using the value in Eq. (|12|) and the ones of light quark 
masses and condensates into the QCD expressions of 



r(4) 

>V,V+A' 



one can deduce the QSSR predictions for 5 



(4) 



and <5y^ given in Table [TJ which agree fairly with the 
range of values from r-decay [14115116] . 
• Four-quark condensate 

Using the ratio in Eq. §§§ and the gluon condensate 
value in Eq. (fT2|) . one can deduce: 

pa s {uu) 2 = (6.4 ± 1.3) 1(T 4 GeV 6 . (15) 

Using this value into the theoretical expression of Sy\ 
one obtains the QSSR predictions in Table Q] (e+e~ col- 
umn), where one can notice a nice agreement between 
the contributions of the D — 6 condensates from the 
LSR in e + e~ and from r-decay Therefore, we can con- 
sider as a final result, their average : 

S ( y ] = (29.4±2.0)1(T 3 => 
pa s {uu) 2 = (4.5 ± 0.3)10~ 4 GeV 6 . 



(16) 



Using the value of a s and (uu) at M T , one obtains: 

p = 3.0 ±0.2, (17) 

confirming the deviations of a factor 2 ~ 3 obtained 
in [7l8l9ll0lllll2ll7ll8ll9l2"0"j from the vacuum satura- 
tion [6] of the value of the four-quark condensate. In 
the V+A channel, the existing results quoted in Table 
UJ are inaccurate except the most recent one from [16[ . 
Considering as the final value, its average with the one 
obtained from Eq. (Til?]) and the expression in Eq. {7}, 
we obtain the value in Table [1] : 

8 ( y ] +A = -(6.1 ± 0.9)10~ 3 . (18) 

• Dimension D — 8 condensates 

There is also a good agreeement in the extraction of 
the D = 8 condensates from the vector component of 
r-decay and QSSR in e + e+ which allows us to take 
the average quoted in Table [1] (e + e~ column). In the 
V+A channel, most of the r-decay results are inaccurate 
but suggest a value much smaller than the one from the 
vector channel. As the theoretical relation in Eq. [7] is 
expected to be (a posteriori) very inaccurate, we shall 

(8) 

consider as a value of <5yj_^, the most accurate estimate 



from [IB] : 

r(8) 
°V+A 



-(0.81 ±0.36)10" 



(19) 



5 A more exhaustive though still incomplete comparison of the 
QSSR estimate of (a s G 2 ) can be found in Table 2 (page 543) of 
the Cambridge book in [7] . In 12.31 , a correlation between the gluon 
condensate and m c has been given. This value would correspond 
to a lighter value of m c which is still in the range of values given by 
PDG 125 1 and from the D-meson sum rules [35 . Further analysis 
to clarify this issue is needed. 



Both results from the V and V±A channels indicate 
large deviations from the vacuum saturation suggested 
by Eq. ([7]). Analogous results have been obtained 
from the LSR and FESR analysis of the V-A channel 
[17118119] . 

• Sum of the SVZ power corrections 

Determinations of the NP corrections from the e + e~ 
sum rules and r-decay data agree, in most case, indicat- 
ing the consistency of the whole picture. However, the 
failure of r-decay for extracting (a s G 2 ), can be mainly 
due to the a 2 suppression of its effect making its extrac- 
tion more delicate. From Table [TJ we deduce: 

8 

5 SVZ ee $ {D) = -(7-8 ± 1.0)10- 3 . (20) 

This result is comparable with the direct fit from 
ALEPH and OPAL [14115] . where the new result is [IB]: 
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£ 

D=4 



5V 



-(5.9 ±1.4)10" 



(21) 



3. Sum of standard power corrections 



We add, to these SVZ power corrections, the ones from 
the pion and ao(980) poles into the longitudinal part of 
the spectral function. It can be written as [36] in units 
of 10~ 3 : 

*r = -16^ 2 4^(l-S|) 2 = -(2-65 ±0.05) 

S a0 = -(0.02 ± 0.01) , (22) 

where we have used f v — 93.28 MeV and f a Q — (1.6 ± 
0.5) MeV [TJ. Therefore, the total sum of the standard 
power corrections reads: 

S st ee S svz + 6^ +6^ + S a0 = -(10.9 ± 1.1)HT 3 . (23) 

4. a s (M T ) from r-decay using the SVZ expansion 

Here, we shall be concerned with the perturbative cor- 
rection which can be expressed in terms of R e + e - 



l + ( j<o) = jL 
MB 



"V'-jayf 1 ^'"-' 24 ' 



The ratio R ee = cr(e + e —> hadrons)/<r(e 



fi+fj, ) 



is expressed through the absorptive part of the correla- 
tor of the electromagnetic current J^: 

Ree = 127iTmn(-i-ie) , (25) 
3Q 2 U(Q 2 ) = i [ d^x e 4 «- a; (0|TJ Al (a;)J^(0)|0), (26) 



with Q 2 = — q 2 . At C(a 4 ), its numerical form for nj 
flavours is : 



R, 



1 + a s + (1.9857 - 0.1152n/) a 2 s 
(-6.63694 - 1.2001371/ - 0.00518^) a 3 
(-156.61 + 18.77 n/ - 0.7974^ 
+0.0215 n)) a 4 , 



(27) 
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where the a 2 ., and af contributions have been com- 
puted respectively in |37l38l39j . The perturbative quan- 
tity 5^ can be evaluated using Fixed Order perturba- 
tion theory (FO) or using the so-called "Contour Im- 
provement" (CI). To order a A s , it reads [1I2I3I4I39] : 



°FO — 



5.202 at + 26.366 ai + 127.079 at 



f (0) _ 



J ci 



= 1.364 a s + 2.54 at + 9.71 at + 64.29 a„ 



(28) 



which, for a reference value a s (M T ) = 0.34, gives : 
4°o = 0.2204 and 5%] = 0.1985 , 



(29) 



where one can notice a slight difference in the evalua- 
tion of 4°] lrom different authors |16l39l4QI41j . This 
difference will only affect the last digit in the esti- 
mate of a s . We use the recent experimental value: 
R ry+A = 3.479 ± 0.011 [US], and the value of the 
sum of standard power corrections 5 s t in Eq. ([23]) . from 
which we deduce: 



5 (0) lst xp - 0.2081 (38) cxp (ll) st . (30) 

Equating Eqs. (J2HD and (JM]) with Eq. (J30D , we obtain 
to order a A s within the SVZ-expansion: 

a s (M T )\ st = 0.3294 (33) cx (10) st (FO) 

= 0.3516 (44) cx (13) st (CI) , (31) 

which, at the Mz scale, becomes: 

a s {M z )\ st = 0.1197 (5)ex(2) st (2)ev (FO) 

= 0.1223 (4) ex (l)st(2)ev (CI) . (32) 

At this stage, the theoretical errors do not take into 
account the estimate of the higher order terms which 
we shall discuss in the next section. One can compare 
these results with some other determinations from r- 
decay within the SVZ expansion given in Tableland the 
runned value at Mz with the one from the Z-width [39] 
and from a global fit of electroweak data at 0(a^) [15] : 



a s (M z )\ N 3 LO = 0.1191 (27) exp (l) th , 
and with the most recent world average 
a s (M z )\ world 

= 0.1189 (10) . 



(33) 



(34) 



One can notice a quite good agreement within the errors 
but, like the ones from [16] and [39], these values are 
on the high side of the world average. Then, one may 
wonders on the possible effects on these results from 
some other (often unwanted) effects beyond the SVZ 
expansion. This will be the subject of the next section. 

5. Power corrections beyond the SVZ expansion 

To the previous standard contributions of the OPE, 
there are also other NP contributions which have been 
overlooked or/and not considered carefully in the exist- 
ing literature. These contributions are expected to be 
present as we work with truncated (non) perturbative 



Table 2 

Different determinations of a s {M T ) from hadronic r-decay to order 
o^, where the difference comes from the variant in the appreciation 
of the non-calculated higher order corrections and of the choice of the 
PT series. ( ) means average of FO and CI and ev means evolution 
from M T to Mz ■ 



PT 



,(Mz) 



R.cf. 



CI 0.3440 (50)ex(70) th 0.1212 (5) eJt (8) th (5) ev \T§\ 

( ) 0.3320 (50)cx(150) th 0.1202 (6)ox(18) t h(3)ov [39] 

FO 0.3156 (30) cx (51) th 0.1180 (4) ox (6) th (2) GV go] 

CI 0.3209 (46) ox (118)th 0.1187 (6) ex (15) ti (3) 8v ED 



This 



FO 0.3294 (33)ox(10) Ht 

CI 0.3516 (44) ex (13) st 

FO 0.3276 (34) ox (86)th 

CI 0.3221 (48)o X (122) th 

( > 0.3249 (29) ox (75) th 



0.1197 (5) ex (2) Bt (2)ev 
0.1223 (4) ox (l) Bt (2) ov 

0.1195 (4) cx (10) th (2) o 
0.1188 (6)e X (15)th(2)o 
0.1192 (4) ox (9) th (2) ov 



Eq. (HTJ 

Eq. gg 

Eq. ([5(3 



Z 

< > 



0.1191 (27) cx (2) th 
0.1189 (10) 
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QCD series which will never describe exactly inclusive 
processes like e.g. r-decay, e + e~ , but instead can only 
give a smearing of it [43] Moreover, even if we are able 
to add the number of terms we want into the (non) per- 
turbative QCD series, we will not be able to get exactly 
the QCD two-point correlator involved in an inclusive 
process such as r-decay and e + e~ because the QCD se- 
ries are factorially divergent. 
• Small size instantons 

Direct instantons are expected to be present in QCD 
for explaining the 77' — it mass shift (the so-called U(\)a 
axial problem [45]). At large Q 2 , it will be highly sup- 
pressed as it can be parametrized by an operator of high- 
dimension D = 9. Its quantitative effect has been dis- 
cussed in previous QSSR literature and has lead to some 
controversy [7j. A phenomenological fit from e + e~ data 
leads to |24"1: 



JV.inst 



20ci v +a ,inst 
-(0.7 ±2.7)10" 



(35) 



Taking for a reference CI perturbative series and the 
corresponding value a s (M T ) in Eq. (|31[) . it would induce 
a shift on the value of a,: 



fo s (M T )| inst = (0.4 ±1.6)10" 



(36) 



which is invisible at Mz within the present accuracy. 
• Quark-hadron duality violation 

A description of the measured spectral function which is 
not possible using PT QCD alone has been initially dis- 
cussed in [44] . This so-called "duality violation" (DV) 



6 As argued in 1441 . lattice calculations may cannot help for solving 
this issue. 
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(deviation of the prediction of the QCD truncated series 
from the measured spectral function) can be quantified 
as gS]: 

ImII(t)|DV = Ke -7 * sin(a + f3t)9{t - i min ) , (37) 

where ft, 7, a and (3 come from the fit of the r-data for 
t m i„ ~ 1.1 GeV 2 <t< M 2 , and depend on the channels 
studied. DV induces an effect [3B]: 

Sdv,v = -(15±9)10" 3 , S dv ,a (2±2)10- 4 , (38) 

respectively in the V and A channels, where 5dv,v 
(which is model dependent) can be relatively large in 
agreement with the rough estimate |<5dv,v| ~ 3% ob- 
tained in [33] 0, while 5dv,a is almost negligible and 
remains to be understood. DV would induce an effect 
in units of 10~ 4 : 

fo s (M T )| DV = 175 ± 101, 8a s (M z )\m = 18 ± 17. (39) 

• Tachyonic gluon from e + e and it sum rules 

This mass can induce a new l/Q 2 -term not present in 
the original SVZ expansion q Some eventual origins of 
this term [33147148149150] and its phenomenological ap- 
plications [52153154117128129130134] have been discussed 
in the literature, as well as its relation with the short 
distance linear part of the QCD heavy quark potential. 
Here, we shall consider that the l/Q 2 -term is purely 
of short distance nature and can mainly emerge from 
the resummation of the infinite terms of the pQCD se- 
ries (UV renormalon). It is expected to provide a phe- 
nomenological parametrization of the unknown higher 
order terms of the PT series as an alternative to the es- 
timates in the existing literature [2|3|14|15|16|55j and to 
the large ^-approximation of the UV renormalon con- 
tribution [33 47 40 . However, its size can depend on 
the order at which the PT series is truncated and may 
(in principle) disappear if infinite terms of the series are 
known [34147] & Its contribution to R T is [53] : 

Mr4 2 iach = ^ 2 4 2 Lh = -2 x 1.05a s A 2 , (40) 

where A 2 is the tachyonic gluon mass estimated from an 
overall fit of the pion sum rule |53|54] and of the 1=1 
part of the e+e~ — > hadrons data [8|52|53] , The average 
of the two determinations gives [54] : 

a s A 2 = -(0.07 ± 0.03) GeV 2 . (41) 

Its presence will only affect slightly the previous deter- 
minations of the condensates having higher dimensions, 
where its correlation to the estimate of (a s G 2 ) from 
e + e~ has been studied in [8]. As originally introduced, 

7 An extension of the model to fit e+e — data above the t mass 
would lead to a slightly lower value of about — (2/3)(6.5±4.2)10 -3 
though consistent with the former within the errors. 
8 An earlier phenomenology of the l/Q 2 -term in r-decay and e+e — 
has been discussed in |5 11521 . 

9 More discussions on its motivation and role in some other 
hadronic channels and QCD phenomena will be reported in a 
forthcoming work. 



it will be, instead, relevant to the estimate of unknown 
higher order PT series, which will be discussed in the 
next section. It would induce, in the V+A channel, an 
effect : 

ChU„o = (46± 20)10- 3 . (42) 

• Tachyonic gluon from large /^-approximation 

It is instructive to compare the previous value of 
the tachyonic gluon (TACH) contribution from the 
one which would be obtained from the large /Jo- 
approximation. In so doing, we take advantage of the 
analysis in [40] . where, for a reference value a s (M T ) = 
0.34, one can predict for the Borel summed PT series: 

<^ 0) = 0.2371 . (43) 

Truncating the PT series known until n = 4, one can 
deduce from Eqs. fl25j! and in units of 10~ 3 : 

Chl^^-^f * 17±0.5(FO) 

~ 39 ± 5 (CI) , (44) 

which agrees nicely with the phenomenological fit in Eq. 
([42]) . We have estimated the error from the deviation of 
the large /^-approximation from the sum of the calcu- 
lated terms of the series truncated at n = 4. However, 
as can been observed from [16140] . though the CI con- 
verges faster than FO, it will never reach the sum of the 
large /^-approximate result, because the CI series sta- 
bilizes for n > 5 — 6 well below the exact value, while 
FO slowly reaches the exact result for n > 7. Therefore, 
the inclusion of TACH with the value in Eq. ([44]) is 
necessary for the PT series to reach the "exact result" . 
In the case of FO, the TACH contribution will decrease 
for increasing numbers of term in the PT series if one 
follows the analysis in 140], while for CI, it will remain 
almost constant. For definiteness, we shall use the large 
/^-approximate result in Eq. (|44|). which can be more 
appropriate and accurate for this specific channel than 
the phenomenological fit. For CI, it would induce an 
effect in units of 10~ 4 : 

<MM T )| tach - ~491±63 , Sa s (M z ) | tach ~ -61±8 .(45) 

• Total contributions 

Adding the previous contributions beyond the SVZ ex- 
pansion, we obtain an estimate of the total "new " con- 
tributions in the V+A channel (in units of 10 -3 ): 

5 nst = 5 inBt + <W + Cch = (2.0±9.4)FO 

= (24.0 ± 10.6) CI , (46) 

where there is a partial cancellation between the DV and 
TACH contributions. These effects can be about a factor 
3 larger than the standard SVZ non-perturbative con- 
tributions given in Table [1] but remain relatively small 
compared to the PT contributions of about 20% which 
we shall discuss in the next section. 
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6. 



5st+nst corrections to a s (M T ) 



Adding the sum of 6 ns t corrections in Eq. (|46[) to the 
value of S s t in Eq. (|23|) . the estimates in Eqs. ([30]) to 
([32]) become to O(af): 



4+nstloxp - 0.2061 (38) ox (ll) st (94) nst (FO) 

= 0.1841 (38) cx (ll) st (94) nst (CI) . (47) 

implying: 

a s (M T ) = 0.3276 (34) cx (10) st (85) nst (FO) 

= 0.3221 (48) ox (14) st (121) nst (CI) . (48) 

However, as explained in the previous section, the result 
from CI is preferred than from FO due to the faster con- 
vergence of the PT series, where the difference between 
the Borel summed series and the known PT series is ex- 

( 2) _____ 

pected to be better reproduced by given in Eq. (|4"2"| 
which remains (almost) constant [4"U] , when higher order 
terms of the series will be available in the future. 

7. Comparison of a s from r-decay and Z-data. 

We run our previous result in Eq. (|48[) at Mz and ob- 
tains: 

a s (M z )\ T = 0.1195 (4) cx (l) st (10) nst (2) ov (FO) 

= 0.1188 (6) cx (2) st (15) nst (2) cv (CI) . (49) 

Both results agree quite well each other and with the one 
from Z-data in Eq. (|33p and from the world average 
in Eq. (13~4"]) . They indicate that the presence of the 
TACH reduces the difference between the FO and CI 
results by a factor 4 [comparison of Eqs. (|3"Tj). ([3"2"j) and 
([49"]) ] , which is not the case of the ones in the existing 
literature. We consider, as a final result, the average of 
the two determinations FO and CI to O (af) and in the 
M 5-scheme: 



(a,(M r )) = 0.3249 (29) cx (9) st (74) nst 
(a s (M z ))\ T = 0.1192 (4) ox (l) st (9) nst (2) 

to which correspond: 

A 3 = 353 (6) cx (2) st (14) nst MeV 

A 4 = 307 (6) ox (2) st (14) nst (2) cv MeV , 

A 5 = 223 (5) ox (l) st (ll) nst (l) cv MeV . 



(50) 



(51) 



In Table [__ we compare our results with the ones ob- 
tained using models which differ in the estimate of the 
strength of the unknown higher order PT terms: 
• The most popular [2I3I14I15I16] is the estimate of the 
coefficient of the next order term either as equal to the 
last term of the PT series or a prediction based on a 
geometric growth of the PT coefficient [2] or using [___] 
the principles of "Fastest Apparent Convergence " (FAC) 
[56] of "Minimal Sensitivity" (PMS) [57] [__ . In [T_] CI 
is preferred over FO as the series converge faster. 



• Using an analogous estimate of the higher order terms, 
Rcf. [35] takes the average of the FO and CI results, 
where 1/3 of the theoretical error (-0.005) comes from 
the estimate of the unknown a b s term. 

• An alternative approach is the use of a toy model for 
the Borel transform of the Adler function beyond the 
n = 4 order. This model favours FO over the CI [40] . 
where the main theoretical error comes from the renor- 
malisation scale dependence of the FO. 

• In [5T], CI has been used in connection with Finite 
Energy Sum Rules having some specific weights. 

• In the present work, all possible (standard <5 st and non- 
standard <5 ns t) QCD power corrections are considered in 
the extraction of a s (M T ). One can notice, from Tabic 
__( a good agreeement between different determinations, 
despite various appreciations on the estimate of higher 
order terms. However, according to the analysis in Ref. 
[46j . none of the previous estimates have carefully con- 
sidered the quark-hadron duality violation (DV) effects 
which cannot be neglected. The DV effects shift the 
results of [16] and [39] to the high-side of the world av- 
erage. In this paper, we point out, that the tachyonic 
gluon mass partially cancels the DV effects and reduces 
the usual systematic difference between the FO and CI 
results by a factor 4 [comparison of Eqs. (|3ip . (|32[) and 

• There is a remarkable agreement between the two de- 
terminations from t [Eqs. ([49]) and ([50]) ] and Z [Eq. 
([3"3"]l ] decays obtained at the same N 3 LO O(a^) level 
of accuracy and at different regions of energy which 
demonstrates the running strong coupling a s of QCD 
predicted by asymptotic freedom. At present, the two 
values of a s (Mz) from r and ^-decays are the most ac- 
curate determinations available today compared to the 
other determinations compiled in 25 42 . Though the 
non-standard power corrections give the largest errors 
of a s (M T ), its accuracy is not largely affected by these 
new contributions. 



8. Tachyonic gluon mass to \V US \ and m s 

• The CKM angle |V US | 

It can be determined accurately [5]: 

|K_ | = 0.2212 (31) oxp (5) th , 

from the observable: 

Rt,S 



|Ks| (Rr/\V ud \i)-5R T ' 
if one uses the experimental values: 



(52) 



(53) 



^ = 3.6661 (12) , — ^ = 0.1686 (47) , (54) 



while the theoretical input only enters into the SU(3) 
breaking quantity: 



SR T = 0.227(54) 



(55) 



'"Some alternative approaches can be found in 



which is one order of magnitude smaller than R T . About 
87% of the error in SR T is mainly due to the value of 
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m s , which has been taken in [53] to be: m s (2 GeV) = 
(96±10) MeV. This error reduces further by a factor 2 if 
one uses the average quoted in Eq. ([5J . One may wonder 
on the TACH contribution to SR T . TACH effects to the 
SU(3) breaking part of the hadronic correlator are [55] : 



mQ 2 ) ( vl A 



sn(Q 2 ) { v ] + A 



l 



16tt 2 



Q 4 

2 fhz 



-f + 16«3) 



8tt 2 Q 4 



(56) 



where C(3) = 1.202056... is the Riemann zeta function. 
Its contribution to 5R T vanishes to leading order in a s 
like other dimension D = 4 operators and can be ne- 
glected in the determination of \V US \ within the present 
accuracy of 5R T . 
• m s from r-decay and e + e~ 

It has been determined from the SU(3)-brcaking r-decay 
widths [5155] : 



R T 



V ud 



\v us 



(57) 



or from r-like decays in e + e using the difference be- 
tween the isovector and </>-meson channels [29|28j : 



Aia = R T ■ 



R T 



(58) 



or with the difference of the usual exponential inverse 
Laplace or/and FESR sum rules [30128] , 
However, TACH contribution is larger in the case of LSR 
and FESR used in |30|28] for determining m s , because 
in these observables, there is no more a cancellation of 
the D = 4 contribution occured in R T . A correlation 
between the values of m s and a s X 2 has been studied. 
FESR analysis with the value of a s X 2 in Eq. ([4"T]) gives 
to order a?: 



• (a s G 2 ) is better determined from QSSR 7 8 24 than 
from r-decay with the robust estimate in Eq. (fT2"]) . 

• In the V channel, r-decay and e + e~ sum rules lead 
to a value of the four-quark condensate pa s (uu) 2 in Eq. 
([16]) which violates the vacuum saturation estimate by a 
factor 3 [Eq. ([17]) ] and confirms previous QSSR results 
from different hadronic channels 7 8 9 1 0|17|18|19|20j . 

• Using the V+A r-decay channel, which is the most 
sensible one for extracting a s , we show in Eqs. (|49[) 
and ([50]) . our estimate to O (af) which takes into ac- 
count the effects of power corrections beyond the stan- 
dard SVZ expansion. Instead of ruining the existing 
estimates without these non standard effects, these new 
contributions improve the agreement of the a s value in 
Eqs. ([49]) and ([50]) with the one in Eq. ([33]) from the 
direct determinations at the Z-mass from the Z-width 
and global fit of electroweak data as well as from the 
recent world average in Eq. Q34[) . The determinations 
of a s (M T ) from r-decay lead to the most accurate value 
of a s (Mz) available today and demonstrates with high 
accuracy the running of the QCD coupling as expected 
from asymptotic freedom ! 

• Finally, we have discussed the role of the tachyonic 
gluon mass by taking the example of the determinations 
of \V US \ from r-decay and of m s from r-decay, e + e~ and 
pseudoscalar channels. Though its effect is negligible in 
the r-like decay sum rules due to the vanishing of the 
D = 4 contribution, to leading order, in this particular 
observable, it plays an important role in the LSR and 
FESR approaches. 
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fh s (2 GeV) 



100 (34) MeV 
93 (30) MeV : 



r— decay , (59) 



in good agreement with the weighted average in Eq. ([5]) , 
but disfavours the value A — to which would corre- 
spond a too small m s mass but having larger errors. 
• m s from the (pseudo) scalar channels 
One should also note that TACH contribution in the 
pseudoscalar pion and kaon sum rule channels reduces 
the estimate obtained for A = to about 5% 1531: 



(m u +rh d )(2 GeV) = 
m s (2 GeV) = 



8.6 (2.1) MeV , 
105 (26) MeV , 



(60) 



which is, in better agreement with the average in Eq. 
([5]). In the same time, it improves the duality between 
the experimental and theoretical side of the sum rules. 

9. Summary and conclusions 

We have reexamined the determinations of the SVZ 
power corrections and their effects on the extraction of 
a s from r-decay data : 
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